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Spin-wave technology (magnonics) has the potential to further reduce the size and energy con¬ 
sumption of information processing devices. In the submicrometer regime (exchange spin waves), 
topological defects such as domain walls may constitute active elements to manipulate spin waves 
and perform logic operations. We predict that spin waves that pass through a domain wall in an ul- 
trathin perpendicular-anisotropy film experience a phase shift that depends on the orientation of the 
domain wall (chirality). The effect, which is absent in bulk materials, originates from the interfacial 
Dzyaloshinskii-Moriya interaction and can be interpreted as a geometric phase. We demonstrate 
analytically and by means of micromagnetic simulations that the phase shift is strong enough to 
switch between constructive and destructive interference. The two chirality states of the domain 
wall may serve as a memory bit or spin-wave switch in magnonic devices. 


Motivated by the aim to reduce energy dissipation in 
electronic devices, spin waves are considered as an alter¬ 
native information carrier in the field of magnon spin- 
tronics [T]. A spin wave acquires a phase shift when it 
passes through a magnetic domain wall (DW) [2]. In 
this Letter, we show that the Dzyaloshinskii-Moriya in¬ 
teraction (DMI) in ultrathin ferromagnetic films makes 
the phase shift dependent on the DW chirality, leading 
to constructive/destructive interference in a two-branch 
interferometer. The mechanism we identify raises the 
prospect of magnonic devices in which DW chirality acts 
as a spin-wave switch. 

It is by now clear that the DMI plays a crucial role 
in the magnetization dynamics of ultrathin films HE], 
due to the broken inversion symmetry at the interfaces. 
The interfacial DMI favors, in a perpendicular-anisotropy 
film, Neel DWs with a fixed chirality [Fig. 0d)] urn , in 
competition with the dipolar interaction, which tends to 
favor Bloch DWs [Fig.[^b)]. The most interesting regime 
is when the two interactions have a comparable strength, 
yielding a DW intermediate between Bloch and Neel Hi 
with two stable minimum-energy configurations (chiral¬ 
ity states) whose in-plane orientations differ by ~ 90°, as 
shown in Fig. [^c). 

Recent experiments demonstrated that DWs can be 
brought into the intermediate regime, and that the DMI 
strength can be fine-tuned by modifying the thicknesses 
of the adjacent nonmagnetic layers i. The internal ori¬ 
entation might be also tuned by an adjacent layer of a 
topological insulator; its surface states induce in the mag¬ 
netic layer an interfacial-DMI-like effect that depends on 
chemical potential and applied electric field [iDHia. 

Our main result is summarized in Fig.j^ where we con¬ 
sider an interferometer in which incoming spin waves are 
divided between two identical waveguides, each contain¬ 
ing a DW. The two DWs are identical in every respect ex¬ 
cept possibly their chirality. When the spin waves rejoin. 
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they are transmitted or reflected depending on the phase 
difference. While it is obvious that spin waves interfere 
constructively if the two DWs have the same chirality 
[Fig-i: a)], we ask if it is possible to achieve destructive 
interference (spin wave blocked) by reversing the chiral¬ 
ity of one DW [Fig. [^b)]. Without DMI, the two chiral¬ 
ity states of the Bloch DW induce identical phase shifts, 
leading to constructive interference [Fig.j^c)]. For strong 
DMI, the DW has a single stable (Neel) configuration and 
the phase shifts are obviously also identical [Fig. [^e)]. 
However, for intermediate DMI, the two equilibrium ori¬ 
entations induce geometric phase shifts differing by as 
much as 180° [Fig. [^d)]. In this regime, the interfer¬ 
ometer can switch between constructive and destructive 
interference - transmission or reflection - depending on 
whether the chiralities are identical or opposite. 

The two chirality states are separated by an energy 
barrier AE (an unfavorable Neel configuration). If AE 
is high enough, spontaneous reversals of chirality due 
to thermal fluctuations are very rare (for the system in 
Fig. we obtain AE = 2.5 x 10“^^ erg = 61 feBTloom)- 
We could consider the intermediate-DMI interferometer 
as a two-state memory device where the transmission 
of spin waves serves as readout mechanism (‘open’ or 
‘closed’). 

Switching does not require modifications of the mate¬ 
rial parameters, nor to insert or remove DWs [2], but only 
to reverse the chirality of one DW, for instance by a field 
pulse normal to the plane of the film. (The ‘field pulse’ 
might alternatively be generated through optomagnetic 
effects dll , provided the light can be focused onto a sin¬ 
gle branch.) The field causes the DW magnetization to 
process as shown in Fig. until, when it is switched off, 
the DW relaxes to the nearest chirality state. 

We have tested the results of Fig. which we de¬ 
rive analytically below, by means of explicit micromag¬ 
netic simulations. The total energy E is given by the 
sum of the usual micromagnetic energy functionals for 
exchange Flex = A //(||i9a;m||^ -I- ||9ym||^) da; dj/, uniax¬ 
ial anisotropy Aani = —K Jf ml da; dy, and dipolar en- 
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FIG. 1. (color online). Effect of the interfacial DMI on the 
magnetization profile m(x) of a DW in a thin film with per¬ 
pendicular anisotropy, (a) Away from the DW, magnetization 
points out of the film (z or — z). Near the DW, the DMI cre¬ 
ates an effective field Hdmi in the — x direction. Depending 
on the competition between the dipolar and DMI interactions, 
the equilibrium configurations OiO *5 and □, shown 

in (b)-(d), are possible, (b) Without DMI, the minimum- 
energy configurations (flux closure) are two equivalent Bloch 
DWs (Oi in dark colors, and 0*i in light colors), whose in¬ 
plane orientations differ by 180°. (c) For intermediate DMI, 
the minimum-energy configurations are intermediate between 
Bloch and Neel. There are two equivalent minimum-energy 
states (★ and ★^), whose in-plane orientations differ by ap¬ 
prox. 90° for an appropriately tuned DMI strength D. (d) For 
strong DMI, a single minimum-energy configuration □ exists: 
a Neel DW with magnetization in the center pointing in the 
—X direction. 


ergy HU, plus a functional 


^^DMi = -2D yy “1 • (Vmz) dx dy (1) 

describing the DMI induced near the interfaces of the ul- 
trathin film . The DMI strength D can be positive or 
negative (we take D > 0 in Figs. and [^. We treat 
the film as effectively two-dimensional (magnetization is 
a function of x and y only), but we do consider the fi¬ 
nite film thickness L in the z direction for the dipolar 
interactions [15] . Here we consider a waveguide made of 
a long strip of ultrathin film with perpendicular magne¬ 
tization {K > 27rMg, where Ms is saturation magnetiza¬ 
tion). The waveguide width W is at least so large that 
the dipolar interactions, in the absence of DMI, favor a 
Bloch DW. 

The interfacial DMI is qualitatively different from a 
DMI oc fff m ■ (V X m) dV present in isotropic bulk 
materials with a chiral crystal structure [IS]. The effect 
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FIG. 2. (color online). Interferometer setup in a thin 
film with perpendicular anisotropy. The two DWs may have 
(a) identical or (b) opposite chiralities. Spin waves enter the 
device from the left. If the chiralities are identical, construc¬ 
tive interference is always obtained on the right-hand side. 
For opposite chiralities, the phase difference depends on the 
DMI strength, as shown in (c)-(e). (c) Without DMI, spin 
waves interfere constructively even if the chiralities are op¬ 
posite ( 0 ) 00 - (d) For intermediate DMI, we find a phase 
difference /S.ip of up to 180° (destructive interference) for op¬ 
posite chirality states (e) For strong DMI, the configu¬ 

rations in both branches are the same (□), trivially resulting 
in constructive interference. In (c)-(e), large arrows repre¬ 
sent the equilibrium magnetization direction m(x). On the 
left, spin-wave basis vectors a, b are defined identically for all 
configurations Oi 0*i D. Their orientation after par¬ 

allel transportation, shown on the right, depends on the DW 
configuration. In (d), notice that the transported basis vec¬ 
tors for ★ and ★' are rotated by 180°. This geometric phase 
difference A:^geom is the dominant contribution to Ay>. 


of a bulk DMI on the interaction of spin waves with DWs 
was considered in Refs. [TB] and [^. Since a bulk DMI 
favors the Bloch DW {Q in Figs. and [^, it does not, 
in the geometry considered here, provide the competition 
with dipolar interactions that is essential to obtain the 
intermediate DW with two equivalent minimum-energy 
orientations differing by approximately 90°. 

Figure shows how spin waves, generated on the 
















3 


y 

I_> 

my 


-0.06 


mx 0.06 


• • ) It u • 


fM'iW/niuin 


tw^niuin 


chirality 

m 


M 






I I I I i f f I« ^ I« 

(a) No DMI 

my 




. • O' 

• • O 

lilH'l I . «t)g| 

I • avg 

constructive 


j I I i ^ ■ ' # • ■ 



t I I I I I • • I I 

(b) Intermediate DMI 


destructive 


FIG. 3. (color online). Micromagnetic simulations of the 
propagation of spin waves in a ferromagnetic thin-film waveg¬ 
uide with perpendicular magnetization (27rM| = 0.907 K, 
L = 31, where I = \J^AJK') through DWs of the indicated 
chiralities I'S- Oi'^)- (a-) Without DMI (D — 0), 

spin waves experience the same phase shift regardless of DW 
chirality, leading to constructive interference (avg = average), 
(b) For intermediate DMI (D = 0.06 A/Z), there is a phase 
difference of almost 180° between spin waves that passed 
through DWs of different chirality, leading to destructive in¬ 
terference. We remark that the attenuation on the right-hand 
side is not the result of Gilbert damping (we take ol = 0.0030), 
but merely represents the present location of the wavefront 


left-hand side of the strip, pass through a DW. We 
solve the Landau-Lifshitz-Gilbert (LLG) equation, for 
relaxed initial states, on a square grid (0.33Z x 0.33Z 
cells) using a self-developed C-l—I- code [TH] with implicit- 
midpoint time integration m- Spin waves are gener¬ 
ated by a space-local, time-periodic in-plane applied field 
(w = 1.70 | 7 |itr/Mg), switched on at t = 0. Each waveg¬ 
uide strip (267Z x lOZ) is simulated in a vacuum-padded 
periodic box (333Z x 271). We calculate the difference 
Aip = (fi' — (fi in phase shift between the two chiralities at 
the right-hand side of the interferometer, comparing in¬ 
termediate DMI to the case without DMI. A phase differ¬ 
ence of up to 180° (destructive interference) is obtained 
for intermediate DMI. 

Fixing AttMs = 3.8 kG and A = 10“® erg/cm [4], 
the other parameters of Fig. become / = 9.9 GHz 
(frequency), W = 127 nm (waveguide width « wave¬ 
length), D = 0.047 erg/cm^, and L = 38 nm, assum¬ 
ing the free-electron gyromagnetic ratio. Spin waves of 
such frequencies and wavelengths can be experimentally 
generated, observed, and visualized Since the 

DMI is an interfacial effect, D is inversely proportional 
to film thickness L [5S]. Extrapolation of the values in 
Ref. [1] (D = 0.5 erg/cm , L = 3 nm) suggests that 
D ~ 0.04 erg/cm^ is realistic for L ^ 38 nm. 

The phase difference Aip between spin waves traveling 


along the two paths (★ and ★') has a geometric [23] 
origin. It is convenient to define (p = (/?geom + </5rei- A spin 
wave causes magnetization to process around its local 
equilibrium direction m(a;) [25] . In the limit of exchange 
spin waves {kx —>■ oo), the dynamics induced by the wave 
is given by the real part of 

m(a;) -b [a(a;) - fb(x)], (2) 

where e > 0 is the infinitesimal amplitude of the spin 
wave (linear regime). The orthonormal basis vectors 
a(a:),b(x) must be perpendicular to m(a:) for all x, so 
that their orientation continually changes across the DW. 
A natural choice is to define a(a::), b(x) according to par¬ 
allel transport, ^ = — (a • ^)m, by which the basis 
vectors, at any given point x, match their orientation in 
an infinitesimal neighborhood of x as closely as possible. 

The function ipy:e\{x) in Eq. determines the phase of 
the spin wave relative to the basis a, b. However, the 
orientation of the basis a, b after parallel transportation 
across the DW strongly depends on the DW configura¬ 
tion (O, O': ★: or □), as shown in Fig. [2])c)-(e). 

This reorientation of a, b implies an additional phase 
shift (pgeom, which is purely geometric in nature. 

It is apparent from Fig. He) -(e) that the geometric 
contribution is approximately given by 


AipgeQjxi ~ 4d, 


( 3 ) 


where i9 is the in-plane angle of the magnetization at 
the DW center, as shown in Fig. [^d). For example, 
we have a geometric phase difference A(/jgeom ~ 180° 
for intermediate DWs with z? = 45°. The value of i? 
is determined by the competition between the DMI and 
the dipolar interaction, as shown in Fig. [^a). While 
in principle A(/3geom depends on the exact shape of the 
equilibrium profile m(a;), we find that the deviation from 
Eq. (§ is at most a few degrees [15]. 

We derive [15] the relative contribution for kx —> 
oo, up to a correction of order as 


D D f°° 


( 4 ) 


where m'y and ruy are the magnetization profiles 'k'y'k 
calculated numerically, taking into account DMI and 
dipolar interactions. Notice that the exchange interac¬ 
tion does not contribute directly to Eq. 0 because the 
basis a(a;),b(a;) (parallel transport) absorbs such a con¬ 
tribution into A(/Jgeom- 

Equation ([^ gives, approximately, 

At^rel = T’rel “ T’rel « ^Wq COS T), (5) 

where wq is a characteristic DW width (wq ~ for 
27rM| ^ K). Notice that At^rei vanishes for D = 0 
(O: O') for large D (□), where z? = tt /2 (Neel wall). 
As shown in Fig. Hb), the contribution of A(/?rei enhances 
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FIG. 4. (color online). Evolution of the magnetization during switching of DW chirality by means of an applied field 
Hz perpendicular to the film in a two-dimensional micromagnetic simulation. Only a part of the waveguide is shown. The 
arrows represent the direction of the in-plane component of magnetization m(f, x, y) and the color the a component. We take 
Hz = 0.09 K/Ms, a = 0.2; other parameters as in Fig. |^b). 



FIG. 5. (a) DW angle i9 as a function of the ratio of DMI 

strength D a nd dipola r interaction, for three film thicknesses 
L, taking ^2-kMq /K = 0.8 |26 |. The angle determines 
whether intermediate DWs (★', ★) are closer to a Bloch or a 
Neel configuration. More “Neel-like” DWs (larger §) are ob¬ 
tained for stronger D. Gonversely, the dipolar interaction pe¬ 
nalizes the Neel configuration (this effect is weaker in thinner 
films), (b) Phase difference Aip between spin waves passing 
through DWs of opposite chiralities (★' vs. ★), as a func¬ 
tion of in the fca, —>■ oo limit, for two values of ^2'KMg/K, 
taking L = 3.01. The dominant contribution Aipgeom is sepa¬ 
rated out. The remaining contribution Aip^ei lowers the value 
of 1 ? needed for destructive interference {Aip = 180°). For 
^2-11 Mq /K = 0.9, Aiprei is larger than for 0.5 because a rela¬ 
tively strong DMI D/ (Awg^) is then needed to obtain a given 
1 ?. 


realized by fine-tuning the DMI strength D [Fig. ^a)]. 

While Eq. ([^ is derived in the short-wavelength limit, 
we have numerically solved the spin-wave normal-mode 
problem [18] for incoming waves of arbitrary wavenum¬ 
ber kx. The phase shifts ip' ^ p depend significantly on kx, 
as in the case without DMI {p = 2arctan (fca,!)”^ |27)). 
but the difference Ap = p' — p between the two chirali¬ 
ties, which is the relevant quantity in our interferometer, 
is weakly wavelength dependent for wavelengths compa¬ 
rable to (or shorter than) the DW width. The weak de¬ 
pendence of Ap on kx can, under certain approximations, 
also be derived analytically [TS]. This observation justi¬ 
fies our approach kx ^ oo. 

In summary, we have shown that the interfacial 
Dzyaloshinskii-Moriya interaction in ultrathin magnetic 
films provides a new way of manipulation of spin waves. 
With this interaction, spin waves experience a differ¬ 
ent phase shift when passing through DWs of different 
chiralities, leading to either constructive or destructive 
interference in a two-branch interferometer. One can 
open or close the transmission of spin waves through 
the device by changing the DW chirality in one of the 
two branches. This opens the possibility of developing a 
memory element or transistor based on the manipulation 
of magnonic currents without charge transport. 


the effect of Atpgeom and merely shifts the critical internal 
angle d for perfect destructive interference {Ap = 180°) 
to a somewhat lower value (more Bloch-like DW). There¬ 
fore the concept of the interferometer spin-wave switch 
is robust: we can always find a DW angle 0° < < 90° 

such that the phase difference Ap between opposite chi¬ 
ralities is exactly 180°. The desired value could then be 
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Supplemental Material 

I. SYMMETRIES AND THE 
DZYALOSHINSKII-MORIYA INTERACTION 

This section provides some background on the distinct 
forms of the Dzyaloshinskii-Moriya interaction (DMI) re¬ 
ferred to in the main text. 


A. Bulk versus interfacial DMI 

For the DMI to be noticeable at the continuum level, 
it is necessary that the central inversion symmetry of 
the system is broken. This may be the case because the 
crystal structure itself is noncentrosymmetric (chiral or 
polar) or because the geometry of the system breaks the 
inversion symmetry, for instance near a surface or inter¬ 
face. A variety of continuum models for the DMI exist, 
which all arise from the same microscopic definition. The 
difference lies in the orientation of the Dzyaloshinskii- 
Moriya vectors Dy that define the interaction between 
atoms i,j. 

We consider two high-symmetry special cases of the 
DMI. The continuum form FIdmi = D f m ■ (V x m) dV 
describes a DMI that is chiral (it changes sign under any 
reflection) but which is otherwise completely isotropic. 
An alternative form Admi = —2i4 / m • V(n • m) dA per¬ 
tains to a system that is polar: the inversion symmetry is 
broken by a director n, which changes sign under central 
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FIG. S6. Directions of the Dzyaloshinskii-Moriya vectors 
Dij and the lattice displacement vectors Vij between atoms 
i,j. (a) For a bulk DMI, the roughly point into the 
direction of Vij. The T)ij “diverge away from” (or “converge 
to”) a given atom i. (b) For an interfacial DMI, the 
“curl around” a given atom i, with a sense determined by the 
normal n. 


reflection; however, the system is achiral in that it re¬ 
mains invariant under any reflections that leave n intact. 
The system is isotropic only in the plane normal to n. 

The interface of a magnetic layer with another material 
defines a polar axis fi (the interface normal), which gives 
rise to a DMI of the polar form (if we neglect any in¬ 
plane anisotropies). For this reason, the polar form is 
usually referred to as the interfacial DMI. The chiral but 
isotropic form of the DMI is usually referred to as the 
bulk DMI, since it is the simplest expression that models 
the DMI in a chiral crystal in bulk (again neglecting any 
anisotropies). A crystal structure is chiral if it cannot 
be superimposed onto its mirror image in any way. We 
remark, for completeness, that a bulk crystal may also 
show an ‘interfadal’ DMI effect if it has a polar crystal 
structure [55]. 

The distinct effects of the bulk DMI and the interfa¬ 
cial DMI are well known from the theory of magnetic 
skyrmions, where it is found that the ‘bulk’ model pre¬ 
dicts Bloch skyrmions while the ‘interfacial’ model pre¬ 
dicts Neel skyrmions [55] . Similarly, it is found that a 
strong interfacial DMI tends to push a domain wall into 
a Neel configuration with a well-defined chirality, while a 
strong bulk DMI stabilizes one particular chirality of the 
Bloch domain wall. 

In Sec. |IB[ we link the cases of bulk DMI and interfa¬ 
cial DMI to atomistic toy models. In Sec. |IC[ we explain 
that the preference of each type of DMI for one or the 
other type of domain wall (Bloch vs. Neel) is a direct 
consequence of its symmetries. 


B. Atomistic and continuum models 


This section shows how the continuum energy func¬ 
tionals for bulk (S9| and interfacial (S12| DMI can be 
derived from the atomistic definition Eq. (S6), given ap¬ 
propriate choices for the Dzyaloshinskii-Moriya vectors 

Di,. 


At the level of individual magnetic moments, the DMI 


is by definition any interaction that can be written as 

Admi = ^ Dy • (m,; X m^), (S6) 


where is the microscopic magnetic moment on site i 
and Dij is the Dzyaloshinskii-Moriya interaction vector 
between sites i and j, which without loss of generality sat¬ 
isfies Dij = —'Dji- For simplicity, let us assume that the 
magnetic atoms of the crystal are arranged as a Bravais 
lattice, where we use the symbol = Vj — Vi to rep¬ 
resent any given (near-neighbor) lattice vector. (Notice 
that, for obvious reasons, the overall crystal structure of 
a noncentrosymmetric crystal must be more complicated 
than just a single Bravais lattice.) 

By translation invariance in the bulk, the magni¬ 
tude and direction of D^- must be a function of 
only. For simplicity, we shall assume that only nearest- 
neighbor interactions are important and that rep¬ 
resents a nearest-neighbor lattice vector. If the ma¬ 
terial is isotropic, the only reasonable choice for the 
Dzyaloshinskii-Moriya vector is = D'iy, where 
D is a (positive or negative) interaction strength. 
The Dzyaloshinskii-Moriya vectors D^- seem to “di¬ 
verge” from any given site i, as shown schematically in 


Fig. S6'a). 

Near a surface or interface, by contrast, the interface 
normal n introduces a preferential direction, which we 
could give a definite sense be defining it to point from ma¬ 
terial A into material B. By symmetry arguments [30], it 
is found that the vectors “curl around” a given site 
i with a well-defined sense induced by the direction of n, 
as shown schematically in Fig. S6T). Considering only 
the highest layer of magnetic atoms below the interface 
and only their nearest-neighbor interactions, the only 
natural choice in the absence of in-plane anisotropy is 
Tiij = Dll X since D^- must be a vector-valued func¬ 
tion linear both in (because of the rule = —Dji) 
and in h (because the polar axis h is the only element 
that breaks inversion symmetry). 


Passing to a continuum theory, Eq. (S61 becomes 


A, 


DMI 


= / [m X (4 • V)m] n 

n. 




(S7) 


where D represents a unit-cell volume and sums over 
all relevant near-neighbor lattice vectors. In the case of 
isotropic bulk DMI with only nearest-neighbor interac¬ 
tions, we substitute D('*.) = Dn.. Equation (S7) becomes, 
in tensor notation, 






(S8) 


where gab = Hn, ^ a^b- Assuming gab = aSab for some 
scalar a > 0 (isotropy), we get 

Admi = -aD J mbebacdarric = -aD J m • (V x m). 

(S9) 
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For the interfacial DMI, on the other hand, we have 
Dij = Dh X where h is the interface normal. For 
simplicity, we take n = z. We get 


^^DMi = ^ J {DeabcZb^c)eademdi^fdfme) n~^cfr 

(SlOa) 

= J DQcfeabc^adeZbmddfme^~^(fr (SlOb) 

= J DgcfiSbdSce - SbeScd)zbmddfmefl~^(fr. 

(SlOc) 


Notice that we now integrate over the interface plane only 
and assume that m does not depend on the perpendicular 
coordinate z (here represents a unit area). Assuming 
Qab = aSab, we get 


£om. = - .,raAm,)nAer (Slla) 


= aD j [(z • m)(V • m) — m • V(z • m)] fl 


(Sllb) 


Using the divergence theorem, we get 
2aD 

lu 


Admi = — - 


n 


[ d^r+^ (f (TO^m)-dzac/, 

Ju “ Jau 


(S12) 

where mz{r) = z • m(r). The boundary term is irrelevant 
if we integrate over all space {U = 


C. Effect on domain-wall energy 

While the bulk DMI is chiral but isotropic, the in¬ 
terfacial DMI is polar but achiral (and isotropic under 
rotations around the z axis). The two forms considered 
here are the only possibilities that satisfy these respec¬ 
tive constraints, regardless of the details of the micro¬ 
scopic model. The same symmetry considerations deter¬ 
mine how the two forms of the DMI respond to Bloch 
or Neel domain walls in a thin film with perpendicular 
anisotropy. In Fig. |S7| we compare the two chirality 
states that exist of either type of wall. The bulk DMI 
selects one of the two chiralities of the Bloch wall as its 
most favorable orientation (energy minimum), and the 
other Bloch chirality as the least favorable orientation 
(energy maximum). On the other hand, the interfacial 
DMI selects one of the two chiralities of the Neel wall as 
its most favorable orientation, and the other Neel chiral¬ 
ity as the least favorable orientation. 

We can interpret this behavior as follows. The two 
chirality states of the Bloch wall are related by a reflec¬ 
tion in the xz plane. Since such a reflection leaves the 
z-director unchanged, it follows that the interfacial DMI 
cannot discriminate between the two chiralities of the 
Bloch wall. The bulk DMI can, because the bulk DMI is 



(b) Neel walls 


FIG. S7. (color online). Symmetries of the two types of do¬ 
main wall, (a) The two distinct states of the Bloch wall are 
related by a reflection that leaves the surface normal n in¬ 
variant. The bulk DMI selects one Bloch state as the most 
favorable and the other Bloch state as the least favorable ori¬ 
entation. The interfacial DMI does not discriminate between 
the two. (b) The two distinct states of the Neel wall are re¬ 
lated by 180° rotation that reverses the surface normal n. The 
interfacial DMI selects one Neel state as the most favorable 
and the other Neel state as the least favorable orientation. 
The bulk DMI does not discriminate between the two. 


not invariant under any spatial reflections. The two chi¬ 
ralities of the Neel wall, on the other hand, are related 
by a 180° rotation around the y axis. Such a rotation is 
a symmetry of the bulk DMI (isotropy) but not of the in¬ 
terfacial DMI (reversal of fi). Asa consequence, only the 
interfacial DMI discriminates between the two chiralities 
of the Neel wall. 

The orientation of a Bloch domain wall will not be af¬ 
fected by a small bulk DMI, as both chiralities represent 
extrema of the bulk DMI. A small interfacial DMI, by 
contrast, causes a change in orientation of the Bloch do¬ 
main wall, reorienting it slightly towards the favorable 
Neel state. 


II. DIPOLAR INTERACTION 

This section specifies the energy functionals used to 
describe the dipolar (magnetostatic) interaction in our 
micromagnetic simulations and numerical calculations. 

Throughout this work, we assume that the magne¬ 
tization is homogeneous in z inside the ferromagnetic 
film (uniform-mode approximation). In other words, 
we assume that the magnetization may be written as 
M(x, 2 /, z) = Msm{x,y, z) = Msll{z/L)m{x,y), where 
L is the film thickness and where n(^) = 1 for |^| < ^ 
and n(^) = 0 otherwise. This simplification can be jus- 
tihed in the regime that L ^ where I is the exchange 
length. 













It can be derived (see, for example. Ref. [3T]) that the 
total dipolar energy (defined per unit thickness) is given. 

in reciprocal space, by 


^ f c\^k 


Edip-2TTMi J TO„(k)5ah(k)m&(k)^2^^2 

(S13) 

with 


guvikx,ky) - {1 Nk) ^2 J 

(S14a) 

9uz{}^x-) ^y) ~ Oj 

(S14b) 


(S14c) 


where k = -\- ky and where the indices u, v repre¬ 

sent the in-plane coordinates x,y. The function ma{r) 
describes the a component of the normalized magnetiza¬ 
tion field. Its Fourier transform rha{k.) is defined as 


ma(k) = / / ma(r)e * d'^r. 


The demagnetizting factor is given by 

1 - 6"'='^ 
kL 


Nk = 


(S15) 


(S16) 


III. PHASE SHIFT IN THE WKB 
APPROXIMATION (fc ^ oo) 


In this section, we provide a derivation of our expres¬ 
sion for the relative part (p^ei of the phase shift in the 
fc —> 00 limit, presented in the main text as Eq. (4). In 
Fig. S 8 we evaluate the geometric phase shift (pgeom and 
the expression (4) for for numerically calculated equi¬ 
librium domain-wall profiles m(a::), confirming the accu¬ 
racy of Eqs. (3) and (5) of the main text. 

An equilibrium magnetization profile m(a;) is found by 
minimization of total energy E under the constraint that 
||m(a;)|| = 1 for all x. It satisfies 


5E 

5va.{x) 


—h{x)ni{x), 


(S19) 


where the expression on the left-hand size denotes a func¬ 
tional derivative. The scalar-valued function h{x) is a 
Lagrange multiplier. The equilibrium profile m(x) suf¬ 
fices to calculate the geometric phase induced by parallel 
transport of the basis vectors. However, to determine the 
relative phase we must solve the normal-mode equation, 
which can be expressed in coordinate-free form [23] as 


We define Nk=o = 1 by conti nuity. Notice that Nk^oo = 
0. We have omitted from Eq. ( |S13| ) some contributions of 
the form JJ||rh(k)|pd^fc = (27r)^ //||m(r)|p d^r, which 
are constant. 

For a system with periodic boundary conditions that 
is defined on a rectangular grid, we evaluate the dipolar 
interaction in O(nlogn) time using a fast Fourier trans¬ 
form (FFT). Finite magnetic elements, such as a waveg¬ 
uides with a finite width, might also, somewhat approxi¬ 
mately, be simulated in this way, provided that sufficient 
zero padding is inserted around the object to minimize 
the influence of the other periodic copies. 

For a one-dimensional magnetization profile m(x), as¬ 
sumed to extend infinitely and uniformly in the y direc¬ 
tion {ky = 0 ), we may write the dipolar energy in real 
space as 

Edip = 27rM| // ma{x')gab{x' - x)mb{x) dx dx' (S17) 


where 


1 / L^\ 

gxx{x) = 5{x) - — log(^l -k , 

(S18a) 

»“W=2iL *“«(' + ?)' 

(S18b) 

gxy{x) = gyy{x) = gxz{x) = gyz{x) = 0 . 

(S18c) 


These expressions can be useful when evaluating the 
dipolar interaction for a system that is not periodic in 
X, such as a domain wall, where we have m(—oo) = z 
while m(oo) = —z. We remark that an efficient numeri¬ 
cal implementation of Eq. (S17) might still evaluate the 
convolution through an FFT, but special care must be 
taken to correctly take into account the asymptotic be¬ 
havior of m(x) for X —)■ —oo and x —>■ oo. 


5'^E 


(5m(x)(5m(x') 


• u(x') dx' -I- h(x)u(x) 


ItI 


iu)[va.{x) X u(x)] 


A(x)m(x), 


(S20) 


where 7 is the gyromagnetic ratio, h(x) is fixed by 
Eq. (S19), and where A(x) is a Lagrange multiplier, un¬ 
der the constraint that m(x) • u(x) = 0 at all x. The 
frequency w is fixed by the wavenumber kx of the spin 
wave far away from the domain wall. 

Equation (|S20|), when written out explicitly in terms 


of u(x) = a(x)a(x) -I- 6 (x)b(x), becomes a complicated 
integro-differential equation. (Nonlocality arises from the 
dipolar interaction.) Here we find an approximate so¬ 
lution using the Wentzel-Kramers-Brillouin (WKB) ap¬ 
proximation, which is exact in the short-wavelength limit. 
We assume that the solution ^ e' ^ locally re¬ 

sembles a plane wave and replace any differential opera¬ 
tor dx with ik'x (analogously for convolution operators) 
in order to solve for k'^ independently for each x. For the 
interactions specified in the main text, including interfa¬ 
cial DMI, we get 


/ Mk'xf + Dmyik'x\ / a\ ^ 

\-^i-Dmyik'x A{k'xf )\k)~ ’ 

(S21) 

where, anticipating the limit \kx \ —>■ 00 , we have written 
only terms that are of at least first order in k'x or kx- 
The exchange interaction acts as a scalar A{k’^Y because 
we define the basis a(x), b(x) according to parallel trans¬ 
port. We substitute the dispersion relation far away from 
the domain wall, a; « ( 2 | 7 |A/Mg)fc^, again up to correc¬ 
tions of constant order. The characteristic equation of 
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FIG. S8. Phase-shift contributions Ai^geom and A^Prei (for 
—>■ oo) obtained from the minimum-energy domain-wall 
profiles, which are numerically calculated as a function of the 
parameters A, K, Ms, L, and D taking the dipolar inter¬ 
action ( |S17[ ) into account. The variables 'd and wo are also 
evaluated, (a) The geometric part Ai^geom follows Eq. (3) of 
the main text (solid line), with a deviation of at most a few de¬ 
grees even in extreme cases (dashed line, ^‘1-rM'^/K — 0.98). 
(b) The relative part A^Srei depends on the DMI strength D 
(as compared to exchange) and on the domain-wall angle i9; 
it follows approximately Eq. (5) of the main text. For given 
Dico jA and the dependence on the third dimensionless pa¬ 
rameter is negligible. The dashed line indicates the relation 
between D and •& if we fix ^2'kMq jK — 0.90 and L = 3.01. 


the matrix in Eq. (S21) has four solutions. We are in¬ 
terested in the solution closest to kx, which is given 

by 


kx{x) = kx + 


1-1 


)• 


(S22) 


The phase induced by the domain wall on top of the phase 
factor is now given in the WKB approximation by 
if = - kx] dx. We find 


V^rel — 


D 


my{x) dx, 


(S23) 


up to a correction of order \kx\~"^. 

Figure [S8] shows a numerical evaluation of the two con¬ 
tributions Atpgeom and Af/jpei to the phase-shift difference 
A(/? = ip' — ip between domain walls of opposite chirality 
{e.g., and ★ in the main text). To describe the behav¬ 
ior of A(/Jrei, we introduce the characteristic domain-wall 


width w = f ■\jTn^ -\- ra^ dx, which we evaluate from the 

equilibrium profile m(a:). The variable wq is then the 
value of w for an equivalent domain wall with D = 0 and 
the other parameters {A, K, Mg, and L) the same. 

For completeness, we derive that the equiva¬ 
lent of Eq. (S23) for a bulk DMI is p^ei = 
-[D/i2A)]j: mx(x) dx. Notice that the latter expres¬ 
sion vanishes for any Bloch domain wall. 


IV. PHASE SHIFT IN A LOCALIZED MODEL 


librium profile and phase shift of a Bloch domain wall. 
In particular, we present a closed-form expression (S41) 
for the phase shift that is valid for arbitrary wavenumber 
kx. The expression suggests that the difference in phase 
shift between the two chiralities is almost independent of 
the wavelength of the spin wave. 

Our analytical treatment complements the WKB ap¬ 
proach of Sec. [m] in two ways. First, it allows one to 
obtain a semianalytical expression for the equilibrium 
domain-wall profile, which in the WKB approach is taken 
as given (ie, calculated numerically). Second, it provides 
an expression for the phase shift for arbitrary kx, where 
the WKB approach considers only the kx ^ oo limit. 
The two approximations made are that we take into ac¬ 
count the main effect of the dipolar interaction as an 
effective local interaction, and that we treat the effect of 
the DMI perturbatively (small D). The results presented 
here are consistent w ith t he WKB expression (Sec. |III[ ) 
if we substitute into (|S23|) the equilibrium profile rn(a:) 


calculated for the approximated dipolar interaction [but 
noti ce th at we need to add the geometric part yigeom to 
Eq. (S23) to obtain the total phase shift]. 


A. Simplified treatment of dipolar interaction 


In our analytical treatment, we follow Ref. |32j in in¬ 
cluding the dipolar effects (shape anisotropy) into a lo¬ 
cal anisotropy energy. We argue that, on scales much 
smaller than the film thickness L, only the ^-function 
part of Eq. ( |S18a ) is im portant and the effect of the 
dipolar interaction dSlTt reduces to a local anisotropy 
J K±mx{x)^ da; with K± = 27rM|. Notice that, even 
though the dipolar interaction is isotropic, the x coor¬ 
dinate plays a special role because we assume that x is 
the domain-wall normal (magnetization is a function of 
X only). The total anisotropy energy is now given by 


^^ani = j{-Kml + Kj_ml) dx. 


(S24) 


where K, K± are positive constants. The in-plane 
anisotropy K± models the dipolar interaction in intro¬ 
ducing a preference for Bloch domain walls (flux closure). 

This localized approximation is formally valid in the 
limit that the magnetization profile m depends only on 
X and extends infinitely not only in the y but also in 
the z direction. In practice, this means that we assume 
normal incidence of spin waves and a wavelength and ex¬ 
change length much shorter than film thickness L. While 
these assumptions are unrealistic in most practical cases, 
they allow us to obtain some analytical results that are 
qualitatively correct. 

In addition to the easy-axis and in-plane anisotropies, 
we take into account the interfacial DMI 

Edmi = —2D / mx{x)m'^{x) dx (S25) 


In this section, we derive analytically, under certain ap¬ 
proximations, the effect of an interfacial DMI on the equi- 


and the usual exchange term Eex = A J'||m'(x)|p dx, 
where a prime denotes a derivative with respect to x. 
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FIG. S9. The function ?(^;s), defined by Eq. ( |S31| |, for 
several values of the parameter s. The special function sech^ 
is shown for comparison. 


B. Equilibrium profile 

In a one-dimensional system, the magnetization profile 
may be described by two functions 0{x), 4>ix), defined by 


file 


6>o(a:) = 2arctan[exp(x/l)], (S29a) 

4>oix) = ±7r/2. (S29b) 


The sign (±) defines the chirality of the domain wall. 
Both chiralities represent equivalent stable energy min¬ 
ima. The positive sign corresponds to the configuration 
as defined in the main text; the negative sign to Q- 
As a result of the competition with Kj_, the interfacial 
DMI (S25) modifies the domain-wall profile, as calculated 
numerically in Ref. [5] . Here, we treat the DMI as a small 
perturbation (|I1| ^ K±l). It can be derived that we get, 
to first order in D, a minimum-energy configuration 


6o{x) = 2 arctan[exp(x/Z)], (S30a) 

d>o(a;) = ±7r/2±, (S30b) 


where s = \/KjKx- Perpendicular magnetization im¬ 


plies s > 1. The function g(^;s), shown in Fig. S9 
uniquely defined as the solution of 


IS 


( rrix 
ruy 


sin Q cos (j) 
sin 9 sin cj) 
cos 6 


(S26) 


The equilibrium magnetization profile of the domain wall 
is a solution of 

= - 2 A 0 " + A(</.'f sin 20 

d9(x) 

+ {K + Kj_ cos^ (/>) sin 29 

- 2D sin^ 9 sin cjjcj)' = 0, (S27a) 

SE 

,,, . = 2A sin 9(j)" + 4A cos 99'(j)' 
d(p{x) 

+ 2K± sin 9 sin </> cos </> 

— 2D sin 9 sin 00' = 0, (S27b) 


where the expressions on the left-hand side are functional 
derivatives of total energy. 

For a system with only exchange A and uniaxial 
anisotropy K [D = K± = 0), it is well known that the 
equilibrium profile, assuming m(— oo) = z and m((X)) = 
—z, is given by 


—cosh^ ^^ sech^ ^ 1 ~ ^ (^31) 


that is even and vanishes at infinity (particular part). 

It is useful to compare Eq. ( S30| ) to the equilibrium 
profile that is obtained if a bulk DMI instead of the in¬ 
terfacial DMI is present. In that case, we obtain the same 
equilibrium profile as in Eq. (S29); in other words, the 


bulk DMI has no effect on the profile of the Bloch domain 
wall. The reason is that, given any Bloch profile, where 
0(x) = ± 7 r /2 and 9{x) is arbitrary, the functional deriva¬ 
tives SEdmi/69{x) and SEdmi/S4>{x) of the bulk DMI 
energy Edmi = D / m • (V x m) dx = —2D f myin'^ dx 
with respect to the profile functions 0 (x), 0 (x) vanish. 


C. Linearized dynamics 


The dynamics is described by the Landau-Lifshitz- 
Gilbert (LLG) equation without damping 


Ms 5m / SE \ 

I 7 I dt ^ dm{x,y)) ’ 


(S32) 


9o{x) = 2 arctan[exp(x//)], (S28) 

where I = A/K is the exchange length, while the func¬ 
tion 00 ( 2 ;) takes an arbitrary constant value. If we set 
ATx > 0, we get the well-known Bloch magnetization pro- 


where Ms is the saturation magnetization, 7 is the gy- 
romagnetic ratio, and E = Eex + Eani + Edmi is the 
total interaction energy. Let us consider small variations 
S9 — 9{t, X, y) — 9o(x) and i50 = 0(t, x, y) — 0o(x) around 
Eq. (S30). We get a linearized equation of motion 


/ _ 7 i (g 2 ^ 52 ) + _ 2sech2 x/l) -[Ms/( 2 | 7 |)] 5 t t (2D/l)E 

[ [Ms/( 2 | 7 |)] 5 tT( 2 D/ 0 Ft -A + d^) + K(1 - 2sech^ x/l) + Kx 


69 

sech{x/l)S(j) 


= 0 


(S33) 
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where 

F = {1/2) sech.{x/l)h{x/l] s)dx cos\i{x/l) 

+ q{x/1] s) ia.-ah.{x/1) (S34) 


We now calculate the effect of the interfacial DMI on 
the transmission phase and amplitude. As above, we take 
into account the DMI to first order in the interaction 
strength D. With some algebraic work, we obtain 


with h = q" — q is the perturbation caused by the DM 
interaction, to first order in D. Notice that F is odd in 
X and breaks reflection symmetry. (The system remains 
invariant under a simultaneous reflection and change in 
polarity, which is equivalent to a rotation around z.) 


( r t' \ _ f 0 e^(vo-vi) 


(S40) 


where ipo is again given by Eq. 


(S39l, and where 


D. Transmission phase shift 


We assume that the solutions are periodic in t and 
y. Equation (S33) becomes a Hamiltonian normal-mode 


problem m- Away from the domain wall (|x| ^ 1), the 
spin waves take the form 


59 

sec\i{x/l)5(j) 


1 

K+Ak^ 

K+Kx+Ak^ 


^i{u]t+k^x+kyy) 


(S35) 


with O’ = (2|7|/Ms )-\/{K + Ak‘^){K F K±^ + Ak'^), 

where k = ■ 

We now calculate the reflection (r) and transmission 
(t) amplitudes of an incoming spin wave with wavenum¬ 
ber kx that propagates in the positive-cc direction. The 
amplitudes r',t' refer to an incoming spin wave wave 
that propagates in the negative-a; direction (wavenum¬ 
ber —kx) and approaches the domain wall from the other 
side. Equation (S35l defines the scattering states. 


For D = 0, the analytic solutions of Eq. (S33) are 
well known [33j . The propagating-wave solutions take 
the form 


]i+Ak^ 1 (S36) 

K+Kx+Ak^ J 

where [27] 

f{x) = —ikxl + tanh(a;//). (S37) 



It follows immediately that 


( 

t' ^ 

1 - 1 

f 0 


\ t 

r' J 

) - 1 

i F'f’o 

0 j 


(S38) 


where 


(fio = arg 


{ikxl - 1 ) 
{ikxl + 1) 


2 arctan - —r. 
kxl 


(S39) 


Notice that the domain wall shows total transmission in 
the localized model. [However, if the full dipolar interac¬ 
tion is taken into account, we find a nonzero reflection for 
long-wavelength spin waves (see also, for example, HZ]).] 


D 


D 


^{K + Ak^){K + K^+Ak^) 

2i^-f + 2Afc2 ■ ^ > 


Notice that only (fi depends on chirality (±). Since (S33) 
defines the scattering problem relative to the basis vec¬ 
tors 6, (j), which have a fixed orientation at x = —oo and 
at X = oo, the expression (S41) includes both the geo¬ 


metric and relative parts of the phase shift. 

Unlike ifo, the chirality-dependent part ifi depends 
only very weakly on wavenumber fc; it is, in fact, almost 
constant in k, as shown in Fig. |S10] We find numerically 
that this conclusion even holds if D is not small or if the 


full dipolar interaction (S17) is taken into account in the 


scattering problem, at least in the regime where the wave¬ 
length is comparable to the domain-wall width or shorter. 
This justifies our approach of taking the kx ^ oo limit 
in the calculation of </3rei in Sec. jlllj (the geometric part 
T’geom is independent of kx by definition). 



FIG. SIO . Dependence of the phase shifts ipo [given by 
Eq. (S39l] and ipi [given by Eq. (S41|[ on wavenumber k. 
Eor a spin wave propagating in the potitive-x direction, the 
total phase shift induced by the domain wall is the sum ipo + 
Notice that i/Pi is almost constant in k. Since only the 
part depends on domain-wall chirality (±), the difference 
in phase shift between the two chiralities - in other words, 
the phase difference obtained on the right-hand side of the 
interferometer shown in Fig. 2(d) of the main text - is also 
almost independent of k. 
































